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ABSTRACT 

Migdal and Kazakov have suggested that lattice QCD with an adjoint repre- 
sentation scalar in the infinite coupling limit could induce QCD. I find an exact 
saddlepoint of this theory for infinite N in the case of a quadratic scalar potential. 
I discuss some aspects of this solution and also show how the continuum D = 1 
matrix model with an arbitrary potential can be reproduced through this approach. 
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1. Introduction 



Recently Kazakov and Migdal have proposed a lattice gauge model 
that, on the one hand, is potentially soluble for large N and on the other 
hand might induce QCD [1]. Their idea was to consider a scalar field, 

in the adjoint representation of SU(N), on a lattice of spacing a, 
coupled to a SU (N) lattice gauge field in the standard fashion, 

S = J2 OTr[-U($(x)) + Y, $(x)iyx)$(x + p)Ut(x)] , (1.1) 

(fi runs over the 2D lattice vectors on a hypercubic lattice of dimension 
D.) Note that there is no kinetic term for the gauge field. This is equiva- 
lent to the infinitely strong coupling limit of the standard lattice action. 
They argued that if one integrates out the scalar mesons (even in the 
case of noninteracting scalars with U($) = |m 2 $ 2 ), then at distances 
large compared to a, one would induce in four dimensions an effective 
gauge interaction, 4^ (In ^^TrF^F^. The basic idea is that the in- 
frared slavery of the scalars, at the size of the lattice spacing, produces 
an effective gauge theory at a larger scale (much larger than the inverse 
scalar mass), which then produces the usual asymptotically free fixed 
point theory. Migdal has analyzed the large N saddlepoint equations 
of this model in great detail [2] and discussed the scaling properties of 
purported solutions, as well as the equations for the spectrum of states 
that emerges [3]. 

There are many problems with this idea. For one the hard gluons 
are not absent and their contribution will overwhelm that of the scalars 
at short distances. Their asymptotic freedom is more powerful than 
the infrared slavery of scalars. It might be possible to overcome this 
problem by considering a general potential U(<&) and fiddling with its 
parameters to arrive at the QCD fixed point. Another issue is that 
the above theory possesses a much larger symmetry than the SU(N) 
gauge symmetry of the usual lattice action. It is not difficult to see 
that, in D dimensions, it is invariant under (D — 1) x [N — 1) extra local 
£/(l)-gauge symmetries. This is because the transformation U^{x)^ 

Vjl(x)U 'y L {x)V ll (x + fia), leaves the action invariant as long as V^{x) is 
a unitary matrix that commutes with Q(x). If V^{x) were independent 
of fi then this would be the ordinary gauge invariance. Thus we have 
D — 1 new gauge symmetries, which are of course isomorphic to the 
special unitary transformations that commute with Thus V^(x) = 
Dn(x)Q(x), where ft(x) is the unitary matrix that diagonalizes $ and 
D^(x) is diagonal. 



This extra symmetry is evident if one considers the naive continuum 
limit of (1.1), by expanding about = 1 + iaA^, since the action will 
not depend on the diagonal components of the vector potential (in a 
basis where $ is diagonal). In fact, one must fix, in addition to the 
usual gauge fixing, the extra (D — 1) x (N — 1) local symmetries to 
eliminate these components. This can easily be done. 

A subset of this symmetry is the local Z^ symmetry, U^(x) — > 

Z^U / j(x)zj 1 , where Z^ is an element of the center of the group. This 
has been noted recently by Kogan et.al. [4]. This symmetry alone pre- 
vents the Wilson loop from acquiring an expectation value, and should 
be broken if we are to recover the QCD fixed point from this formula- 
tion. Presumably the same is true of the extra continuous symmetries 
discussed above. 

In this paper I shall construct an exact solution of the saddlepoint 
equations for the case of a quadratic potential. Remarkably, it turns out 
that the eigenvalue distribution, for large N, is given by a semi-circular 
law with a mass parameter that depends on m and on the dimension of 
space-time D. This trivial solution bears no resemblence to continuum 
QCD. Instead it describes the strong coupling lattice theory- with color 
singlet mesons. 



2. The Quadratic Potential 

The model described by (1.1) can be reduced to an integral over N 
degrees of freedom per site, by integrating out the matrices that diag- 
onalize $(x). This can be done, link by link, in terms of the Itzykson- 
Zuber integral (A($) = Ui<ji^i " [5], 

'(•^-/w/e^I^-c^g. (2.1) 
The effective action for the eigenvalues is then, 

S = Y}- N U($i(x))]+52 hi\$i(x)-$j(x)\^HI($(x)Mx+w)l 

x,i x,i^j x ^ 

(2.2) 

In the large N limit the integral will be dominated by a translationally 
invariant saddlepoint for the density of eigenvalues, p{y), of the matrices 
The saddlepoint equation is 



p 



J dv J^L = - DL i mjv _ [WH] | M . (2.3) 



Migdal has simplified these using the Schwinger-Dyson equations that 
are satisfied by [2]. These are consequences of the fact that 

/ satisfies Tr[(^) k ]I = Tr(*) k I. The net result is that one derives 

an equation for the function V'{z) = f dz^-, whose imaginary part is 

lmV'{v) = -irp(v), 



Rey'(A) = P I ± H *--^'M-f ReV' M + t Mv) 
' 2711 ^--^T)U'{u)-^ReV'{u)-nip{uy 



Migdal has studied the scaling properties of conjectured solutions 
of this equation. Here we shall find an exact solution in the case of a 

quadratic, Gaussian, potential, U(&) = |m 2 $ 2 . First, we note that this 
equation simplifies dramatically for D = 1. This is not surprising since 
in one dimension the gauge field plays no role, can be reabsorbed into 
a definition of and (1.1) reduces to the standard action for a scalar 
field in one dimension. In particular for the quadratic potential the path 
integral is a Gaussian, 

Z = J n^-^^tf^^ 11 • (2.5) 



Thus the eigenvalues of $ will be given by the semi-circular distribution, 

/ 2 2 

namely np{v) = y fi — ^j-, where fi is determined by the mean of the 

squares of the eigenvalues, (^Tr($ 2 )) = ^. It is therefore sufficient to 

calculate the expectation value of -^Tr(<I> 2 ), which is given by the one 
loop integral, 



I Tr (^) = / * * = - 1 = I . (2.6) 

N ' J 27rm 2 + 2cosp Vm 4 - 4 M 

— 7T 



It is easy to verify that this solves (2.4), using the fact that 



V'{z) = ^-^^--^ ReV'(u) = >, for \v\ < A . ( 2 .7) 

However, if we return to (2.4), we see that the integral involved is of the 
same form for any D, as long as ~ReV'(v) is linear in v. This suggests 
that we can find a solution of (2.4) with a semi-circular distribution of 
eigenvalues for a quadratic potential in any dimension. 

We use the result that, 



2 



(2.8) 

(z—-) 

This result is easily established (change variables v = — i ^= , inte- 
grate once by parts and then the z integral can be done by contour 
integration). In our case, assuming the semi-circular distribution in 

(2.7), b = m 2 + ^ D D ^ fi. On the other hand, from the integral it must 

be that b 2 = fj 2 + 4. Therefore we derive a quadratic equation for /i, 
whose solutions are, 



m 2 (D-l)± £»y/m2-4(2D-l) 
M±(D) = ' (2 ' 9) 

In particular for D = 1, = = Vm 4 — 1, which agrees with the 

direct calculation of (2.6). 

Thus we have found a saddlepoint for any D, as long as the potential 
is quadratic. It is given by (2.7), with ji given above. If we start with 
a large mass, m 2 > m 2 = 2D, then only fi + is positive and must be 
chosen. This solution can be verified to satisfy the large mass expansion 
discused by Kazakov and Migdal [1]. I have explicitly checked that a 

semi-circular distribution of eigenvalues, with \i = = m 2 — ^ — 
^§nr— - — 4P ^io~ 1 ' > • • •, satisfies the expansion up to terms of 

fit fft fit 

order (^) 14 . 

The free energy, F = — N 2y Q i m <Z\ can be calculated at the saddle- 
point by the following device. We use the fact that for D = 1 the free 



energy can be evaluated by Gaussian integration to yield F(D = 1) = 

gln[m 2 + 2cosj9] = Jln(^±^S). Then we note that for any 
D the free energy is given at the large N saddelpoint by, 

F(D) = ^JlY<E> s 2 - |^A 2 ($ S ) - plnl($ s ,$ s ) , (2.10) 

where $ s is the saddlepoint scalar field. The first two terms are easily 
evaluated for a semi-circular distribution of eigenvalues, which allows 
us to express the logarithm of the Itzykson-Zuber integral, for a semi- 
circular distribution of eigenvalues with parameter /i, as 

^ln/O) = v/ ^+ i ~^ - ^ln( v ^y +/i -). It then follows that the free 

energy, for a semi-circular distribution of eigenvalues with paramter \i 
is 

It is easily checked that the saddlepoint of (2.11) determines \i to satisfy 
(2.9). 

For large mass it is easily verifed that F is bounded from below, 
and achieves its minimum at \i = //+. For D < 1, //+ vanishes at 

m 2 = 2D, and thus we can take a continuum limit of the theory. For 
D > 1, ji+ is always non- vanishing. Instead as m 2 — > 2D the other 
branch vanishes. However, we must approach the critical value of 

m 2 from below, since only then is /i- positive. If we examine the free 
energy in this case we find that it is unbounded from below. Indeed 

F(D) m ~ 2D . The saddlepoint at \i = //_, which vanishes as 

m 2 — > 2D — 0, is a local maximum. Thus although we may construct a 
continuum limit about this point it might contain tachyonic excitations. 

It is not surprising that F is unbounded from below, for m 2 < 
2D, since the action can be bounded, for translationally invariant $, 

by S = iVVol.Tr[-^$ 2 + ^ QU^UU > NVol.Tr[-2i-V + D$$], 

which is unbounded from below. Thus, this is not a stable region of the 
parameters. Unless there exist other, healthier, large TV saddlepoints of 
the action, then there there would not appear to be a scaling solution 
that would give a continuum theory for a Gaussian potential for D > 1. 
Nonetheless, the saddlepoint we have found will describe the large N, 
strong coupling limit of QCD with adjoint matter. 



3. Arbitrary Potential in One Dimension 



Now let us consider the one dimensional system for an arbitrary 
potential, namely the model described by 



This model is hard to solve, except in the double scaling limit, where 
it was solved exactly as long as the spacing between the points is not 
too big [6]. However, in the continuum limit, the model can be solved 

in terms of N free fermions in the potential ^m 2 $ 2 + V(<&), and the 
density of eigenvalues is easily calculable in the large N limit. To take 
the continuum limit we must scale, 

$ n ^^; m 2 = 2 + a 2 M 2 ; [/($)= aW{ip) . (3.2) 

Oj 



Then the above partition function reduces, as the lattice spacing a — > 0, 
to the continuum integral, 

Z= / > p^(^) e -^/^[ Tr ^^) 2+iy (^))+i( a ^( t )) 2 ] . (3.3) 



Now, this model is easily solved, following [7], by reducing it to the 
problem of N free fermions, (the eigenvalues of <£>), with the one body 

2 2 

Hamiltonian, H = ^- + ^|-A 2 +M / (A), with playing the role of Planck's 
constant. In the WKB approximation [i.e. N — > oo limit), 

P{E} = nIe = J ^r s[E ~ t ~ ~ x ~ w(x)] 

I M 2 

-> p(x) = \/2(e - —x 2 - iy(x)) . 



The Fermi level, is determined by the condition that J^ x dxp(x) = 1. 

Let us return to the matrix chain, and write the equations for this 
case, 



A-(l + ^-jai^) + W ] 
Imy» = -7170(1/) 

(3.5) 

Now rescale the variables so that A = is = p(y) = y / a/(z>), 
y'(A) = aF(X). In terms of these the equations read, 



Re-F(A) = 1 / ^ !n[ A - * ~ + , , 

a J 2m v \-y-oW i> _^ 2 W '{p)-iaT,f{vy . (3.6) 

ImF(z>) = -tt/(z>) 



Expanding (3.6) in powers of a we find 

ReFCx)= ( d ,m + i[j M ^ +w '^ 



X — P 2 J (A - s>) 2 

(A - P)' 



— - / eh/ ' \; + 0(a 3 ) 



(3.7) 



Now as a — > we see that these equations are trivially satisfied for 
any function f{v). We now demand that the term of order a 2 vanish. 
Integrating by parts the last term this condition can be written as 

T / ([iW 2 £> + V^ / (z>)]/(f>) + 7T 2 / 2 (£>)/ / (z>)) . Consequently M 2 z> + 

W f (i>) + 7T 2 f (i>) f f (v) = 0, an equation that we can solve for /, 



tt/(z>) = ^2E - M 2 v 2 - 2W{y) . (3.8) 



E is a constant that can be fixed from the normalization of p. This 
agrees with the fermionic solution (3.4). 



4. Conclusions 



We have found that infinite coupling QCD, with scalar adjoint mat- 
ter, is dominated for large N by a very simple saddlepoint. The eigen- 
values of the scalar field are distributed according to the semi-circular 
law. This is quite remarkable, especially as the effective scalar action is 
certainly not a Gaussian. It will be very interesting to evaluate the fluc- 
tuations about his saddlepoint-ie. to work out the effective Lagrangian 
describing the confined mesons in this lattice gauge theory. We note 
that for zero gauge coupling the eigenvalues also have a semi-circular 
distribution. In this case \i will be given by the analog, in dimension D, 
of (2.6). The fact that the eigenvalue distribution is semi-circular for 
both zero and infinite coupling makes one wonder whether this could be 
the case for finite coupling as well! 

What have we learned about induced QCD? We have certainly 
learned that the program cannot work for a Gaussian potential. This 
is consistent with the analysis of the scaling solutions of the genral 
equations [8]. So if the idea is to work one must consider non-trivial 
potentials and hope that the non- asymptotically free scalar interactions 
do indeed induce QCD. Although I do not see why this could not happen 
the physical mechanism is mysterious. 
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